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ABSTRACT 

We discover a faint shell-type radio and X-ray source, G353.6-0.7, associated with HESS J1731-347. 
G353.6-0.7 is likely an old supernova remnant (SNR), based on radio (0.8 GHz, 1.4 GHz and 5 GHz), 
infrared (8 /zm from the GLIMPSE Legacy Project and 21 /zm from the Midcourse Space Experiment), 
and X-ray (0.1 keV - 2.4 keV from the ROSAT survey and 5 - 20 keV from the INTEGRAL survey) 
data. The SNR, centered at (7, b)=(353.55, -0.65) with a radius of ^ 0.25°, closely matches the outline 
of the recently discovered extended TeV source HESS J1731-347, which has no previously identified 
counterpart. A diffuse X-ray enhancement detected in the ROSAT all-sky survey is coincident with 
lower half shell of the SNR. Therefore the SNR is the best radio counterpart of both the HESS source 
and the diffuse X-ray enhancement. G353.6-0.7 has an age of ^ 27000 yrs. Altogether, the new 
discovery provides the best case that an old SNR emits TeV 7-rays. : 

Subject headings: (ISM:) supernova remnants-X-rays :-gamma rays: observations-radio contimmm: 
galaxies-radio lines: galaxies 
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1. INTRODUCTION 

About 52 Galactic TeV 7-ray sources have been de- 
tected by ground-based high energy telescopes (e.g. the 
High Energy Stereoscopic System (HESS), Aharonian et 
al. 2006; the MAGIC telescopes, Albert et al. 2006) 
so far. Some of them do not have counterparts in other 
wavebands. Seeking their Galactic counterparts in other 
wavebands could play a key role for us to understand the 
origin of these sources and distinguish different acceler- 
ation mechanisms to emit TeV 7-rays. Among identi- 
fied counterparts, young shell-type supernova remnants 
(SNRs) are one of the four Galactic populations to gen- 
erate very high energy (VHE) 7-rays, i.e. pulsar wind 
nebulae (PWN), X-ray binaries, shell SNRs and a young 
stellar cloud Westerlund 2 (Aharonian et al. 2007, Landi 
et al. 2006, Torres et al. 2003, Enomoto et al. 2002). 
Recently, the old SNR W41 has been proposed to as- 
sociate with the extended TeV source HESS J1834-087 
based on observations (Tian et al. 2007; Leahy & Tian 
2008), which is theoretically expected (Yamazaki et al. 
2006, Fang & Zhang 2008). However, this is the only 
likely association of an old SNR with a TeV source up 
to now. We report here the discovery of a new faint ra- 
dio and X-ray shell which closely matches the extended 
TeV source HESS J 173 1-347 without previously identi- 
fied counterpart (Aharonian et al. 2008). We use new 
radio, infrared and X-ray data to identify the shell as an 
old SNR, therefore obtaining second association of a TeV 
source with an old SNR. 

2. RADIO AND X-RAY OBSERVATIONS 
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The 1420 MHz radio continuum and 21 cm HI emission 
data sets come from the Southern Galactic Plane Survey 
(SGPS) which is described in detail by Haverkorn et al. 
(2006) and McClure-Griffiths et al. (2005). The SGPS 
images the HI line emission and 1420 MHz continuum 
in the range 253° < I < 357° and |6| < 1.5°, mostly in 
the fourth quadrant of our Galaxy, with the Australia 
Telescope Compact Array (ATCA) and the Parkes 64m 
single dish telescope. The continuum observations have 
a resolution of 100 arcsec and a sensitivity better than 

1 mjy/beam. The HI data have an angular resolution of 

2 arcmin, a rms sensitivity of ~1 K and a velocity reso- 
lution of 1 km s~^. We combined the ATCA and Parkes 
data to estimate the flux density of the SNR, because 
the Parkes data properly includes scales of around 0.5° 
as well, where the ATCA becomes less sensitive. 

The vicinity of SNR G353.6-0.7 was observed by the 
Position Sensitive Proportional Counter (PSPC) on the 
ROSAT X-ray Observatory during the all-sky survey in 
1993 for a total source integration time of 300 seconds 
(RS932341N00). We obtained an exposure-corrected in- 
tensity image of the remnant in the 0.1-2.4 keV band, 
smoothed with adaptive Gaussian kernel to achieve 
Signal-to-noise ratio ^ 5 (see Fig. 1). 

3. RESULTS 

3.1. Continuum Emission 

The new SNR has a nearly-circular shell-type structure 
with angular size of 0.5° centred at (i, b)=(353.55, -0.65) 
in radio continuum (see the first row of Fig. 1). A bright 
compact HII region G353.42-0.37 is adjacent to the west 
of the SNR. The 1420 MHz image reveals a clear fila- 
ment along its south boundary. The 843 MHz image from 
the Molonglo Observatory Synthesis Telescope (MOST, 
Green et al. 1999) also shows the filament. The filament 
of the SNR is also seen at 5 GHz by Parkes observations 
(Haynes et al. 1978), but not as clearly asat 1420 MHz or 
843 MHz due to the lower spatial resolution (4.iP"™e). 
A few compact sources within the remnant are resolved 
in both SGPS and MOST maps. Specifically, the bright 
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source G353.45-0.68 is close to the center of the SNR. 
The SGPS image has a lower resolution (100 arcsec) but 
higher sensitivity (<1 mJy/beam) than the 843 MHz 
map (43 arcsec and 1 to 2 mJy/beam). The SGPS image 
shows more faint emission than the MOST map because 
the MOST image misses short-spacings more severely 
than for the SGPS. The 1420 MHz and 843 MHz images 
have similar shell-type outlines for the SNR. 

21 cm weak polarized emission is seen toward the south 
filament of SNR G353.6-0.7. However, smoothing the 
same data with a resolution of 5 arcmin to enhance any 
extended features shows no convincing polarized emis- 
sion coming from SNR G353.6-0.7. Absence of observ- 
able polarization is consistent with the expectation of 
complete depolarization at a distance of >~1 kpc toward 
the Galactic centre. 

We derive an integrated flux densitv of 2.2±0.9 Jy for 
G353.6-0.7 at 1420 MHz by a new' method. We az- 
imuthally averaged the radio intensity in rings around 
the center of the remnant. Only one half of the rem- 
nant, the part at lowest Galactic latitudes, was used be- 
cause the half closest to the Galactic plane is confused 
with other sources. The integrated flux density has ex- 
cluded the contribution from a few compact sources (e.g. 
G353.45-0.68) within the remnant. Fig. 1 (lower right) 
shows radio flux as a function of radius of the remnant in 
diamonds. The background was estimated as the aver- 
age intensity at radius R > 0.4° and subtracted from the 
averaged flux. The profile of the remnant was estimated 
by a polynomial. This profile was integrated over radius 
from the center of the remnant out to i? = 0.37° and 
over 27r in angle, which resulted in a total flux density 
estimate of 2.2 ± 0.9 Jy. The 843 MHz data misses too 
many low spatial frequency components to be used to ob- 
tain a reliable flux density. However, we obtain an inte- 
grated flux density of 1.4±0.3 Jy at 5 GHz by same above 
method, giving a spectral index of a=0.4(^5 (S^ ~ 
consistent with a non-thermal feature. 

3.2. HI Spectra and Distance Constraints to 
G353.6-0.7 

The compact HIT region G353.42-0.37 is adjacent to 
the west of the SNR. Fig. 2 shows the HI absorption 
spectrum features in the range of 9 to -24 km/s which 
is derived using our improved methods (Tian, Leahy & 
Wang 2007; Tian & Leahy 2008; Leahy & Tian 2008). 
Lack of absorption for velocities less than -30 km/s shows 
G353.42-0.37 is at the near distance (3.2±0.8 kpc, using 
a flat Galac;tic rotation curve model and the most re- 
cent estimates of the parameters for this, i.e. velocity 
V/{=Vo=200 km/s and Ro=7.6 kpc, see Eisenhauer et 
al. 2005) of its recombination line velocity (~ -16 km/s). 
Based on the fact that a number of HESS sources are 
detected where high-energy particle accelerators are lo- 
cated near H H regions, and other such accelerators with- 
out proximate sources of low-energy photons are weak or 
undetected at TeV energies, Helfand et al. (2007) sug- 
gested that TeV emissions of HESS sources could orig- 
inate from inverse Compton scattering of the; starlight 
from the nearby H II regions. Given the proximity of 
bright HII regions to several HESS sources, SNR G353.6- 
0.7 may be at ^ 3.2 kpc, same as G353.42-0.37. We have 
also obtained the absorption spectrum of compact source 
G353.45-0.68, and derived that this source is likely an ex- 



tragalactic source due to absorption features in the whole 
negative velocity range. Therefore, we exclude a possibil- 
ity that G353.45-0.68 is a pulsar or pulsar wind nebula. 

3.3. Infrared and X-ray Emission 

We analyze new high resolution data (images) at 8/im 
from the Galactic Legacy Infrared Midplane Survey Ex- 
traordinaire (GLIMPSE) Legacy Project (Benjamin et 
al. 2003), and a mid-IR image at 21 /xm from the Mid- 
course Space Experiment (Price et al. 2001). These im- 
ages (all are centred at [b=353.6, l=-0.7]), show no ther- 
mal emission down to a limit of 0.4 mJy (5 a) at 8/im and 
200 mJy at 21/im associated with G353.6-0.7. This sup- 
ports above-mentioned conclusion that G353.6-0.7 has a 
non-thermal emission character and is a new SNR. 

The lower left panel of Fig. 1 shows the ROSAT X- 
ray image of SNR G353.6-0.7 in the the soft band (0.1-2.4 
keV). The diffuse X-ray emission appears only within the 
lower-latitude half of the radio remnant. The coincidence 
of the X-ray morphology with the radio morphology (two 
bright X-ray sources overlap the radio filaments along 
SNR-shell, and both the X-ray and the radio emission 
have same shape) suggests that they are physically asso- 
ciated. The fact that X-rays arc detected from only the 
lower half is likely due to absorption blocking X-rays from 
the upper half. This is supported by the available low 
resolution 12C0 maps of this area (Dame et al. 2001), 
showing increasing column density toward the Galactic 
plane. 

The estimated unabsorbed X-ray flux of SNR G353.6- 
0.7 is ~ 8 X 10~^^ ergs cm~^ s~^ in the energy band of 
1-2 keV (the energy band below 1 keV is affected by 
absorption) by fitting the ROSAT spectrum with a ther- 
mal plasma model (i.e. XSPEC MEKAL with a fitted 
temperature kT > 2.8 keV). This flux is corresponding 
to a luminosity of 9.8x 10^'^ ergs s~^ assuming a dis- 
tance of ~ 3.2 kpc. Due to low number of X-ray counts, 
this model is quite uncertain. The ROSAT spectrum is 
extracted from the lower half of SNR G353.6-0.7 (half 
circle with radius of 15') centered on the remnant. 

4. DISGUSSION AND CONCLUSION 

4.1. An old SNR 

SNR G353.6-0.7 has a large angular size and a large 
physical size (^ 28 pc in diameter for a likely distance 
of ^ 3.2 kpc), and faint radio and X-ray emissions. This 
suggests that G353.6-0.7 is quite old, because young 
SNRs (< 2000 yrs) usually have a size < 10 pc, and 
produce strong radio and X-ray emission. Applying the 
Sedov model (Cox 1972), we estimate its age based on 
the known radius R (14 pc) and ambient density no as- 
suming a typical explosion energy (10"''^ tirg). The likely 
absorption of X-rays from the western half of G353.6- 
0.7 hints that ambient density no is not very low. The 
density no can be estimated from the 7 flux of HESS 
J1731-347 (F^«l. 2x10-1^ cm-^ s"^ which is given by 

dF^/dE{> Ermn) = N„ {tTcvT^ ' Here Emin=Q.^ TeV, 
r=2.26±0.10, A^o=6.1±0.8 xlO^^^ ^.^^^-2 g-i g^,^ p^^^a- 
ronian et al. 2008). Using the formulae 1 and 2 of Aha- 
ronian et al. (2006), we obtain the no of ~ 5 cm~^ (as- 
suming a conversion efficiency of 0.1 - 0.2 and explosion 
energy 10^^ erg). If G353.6-0.7 is still in the Sedov- Taylor 
phase, it has a Sedov age of ~ 25000 yrs. However, the 
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Fig. 1. The first row shows the continuum images of G353.6-0.7 at 1420 MHz (left) from the SGPS (ATCA data only) with a beamsize 
of 100 arcscc and at 843 MHz (right) from the Molonglo Galactic Plane Survey with a bcainsizc of 43 arcscc. The left of the second row 
presents the ROSAT/PSPC image. The right is azimuthally averaged flux at 1420 MHz over the low-Galactic-latitude half of G353.6-0.7 as 
a function of radius from the center of the remnant at [5=353.55", l=-0.65'']. The horizontal solid line denotes the background, estimated 
from the intensity at i? > 0.4°. The wave solid line is a polynomial fit to the radial profile of the remnant. The direction of North (N) and 
South (S) is marked on the first plot. 
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Fig. 2.— The HI spectra of HH G353.42-0.37 



radius at which cooling affects the dynamics of an SNR 
is ~ 10 pc for no ~ 5 cm~^, so G353.6-0.7 has already en- 
tered the radiative phase. So wc here obtain a radiative 
age of - 27000 yrs for G353.6-0.7. 

4.2. Radio, X-rays and TeV j -rays from G353.6-0.7 



The observations from radio, X-ray and 7-ray all in- 
dicate that SNR G353.6-0.7 is the counterpart of the 
newly detected TeV 7-ray source HESS J1731-347 (see 
Fig. 3). The TeV 7-rays overlap the entire radio SNR, 
including the east half seen in X-rays and the west half 
not seen in X-rays. The ^^CO column density increases 
toward the Galactic plane, consistent with absorbing 
the X-rays from the west half of the SNR, but allow- 
ing both radio and TeV 7-rays to be observed. Old 
SNRs emitting TeV 7-rays are rare. W41 has been 
considered as the only likely candidate so far (Tian et 
al. 2007). Young SNRs found to show 7-ray emis- 
sion have F^:^i_ioTeV)/Fx{2-wkeV) of typically below 
~2, -FIy(i_ioTeV)/-P'ra*o(o.oi-iooGHz) of typically below 
~10 (Yamazaki et al. 2006). 

The flux density of G353.6-0.7 is Si42aMHz~ 2.2 
Jy, giving an integrated flux from 10^ to 10^^ Hz of 
Frodio=5.2 X 10~^^ ergs cm~^ s~^, assuming a radio 
spectral index of ~ 0.5. We calculate the 7-ray flux of 
F7(i-ioTey) ~ l-7x 10"^-^ ergs cm"^ s^^ in the 1-10 
TeV band, so the ratio R = F^/Fradio is ~ 33. Also 
given a distance of 3.2 kpc to G353.6-0.7, the 7-ray lu- 
minosity is ~ 2x10^^ ergs s~^ in the 1-10 TeV range, 
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Fig. 3.— The 1420 MHz image of G353.6-0.7 with 
ROSAT image contours (red: 0.11, 0.19, 0.23, 0.8 in units 
of 10~^countss~^pixel~^) and HESS J1731-347 image contours 
(green) overlaid. The HESS J1731-347 image has same contours as 
Fig. 5 of Aharonian et al. (2008), i.e. starting at 4a in la steps, 
which is of 7-ray excess counts smoothed with a Gaussian filter 
with standard deviation 0.1° 

which is larger than that of RX J1713. 7-3946, one of the 
intrinsically brightest galactic TeV Gamma-ray sources 
discovered and thought to originate from pion decays. 
So, SNR G353.6-0.7 is a powerful cosmic ray accelerator. 

The theoretically predicted 7-rays to X-rays flux ra- 
tio for old SNRs is possibly >'-~^ 100 (Yamazaki et al. 
2006), giving an upper limit to 2-10 keV X-ray emission 
from G353.6-0.7 of about 1.7x 10"" ergs cm'^ 3"^ We 
checked the INTEGRAL survey which covers G353.6-0.7 
and found no detection within the SNR (JEM-X in en- 
ergy range of 5 - 20 keV). Assuming an individual IN- 
TEGRAL observation for an integration time of 2 ks, the 
5 a source detection limit is 1.2 x 10"^" ergs cm~^ s~^ 
on-axis in the energy range of 2 - 10 keV (Brandt ct al. 
2003). This is consistent with the derived X-ray upper 
Hmit flux from the old SNR. 

The fact that the X-rays are visible in SNR G353.6- 
0.7 shows that SNR shock velocity is still fast. This 
means that G353.6-0.7 is not too old (i.e. <10^ yrs), be- 
cause an old SNR is expected to emit few thermal X-rays 



which originate from the high temperature gas heated by 
the fast-enough shock. Because of poor statistics in the 
ROSAT X-ray spectrum, we can not exclude that part 
of the X-rays are from synchrotron emission. Because no 
pulsar/PWN is detected within the remnant until now by 
radio and X-ray observations, it is unlikely that the TeV 
7-rays originate from particles accelerated to high energy 
at a termination shock formed in the pulsar wind. This 
strengthens our conclusion that the TeV 7-rays originate 
from the old SNR. 

TeV 7-rays can originate from an old SNR or from a 
molecular cloud encountered by the SNR, via pion de- 
cay with the pions produced by proton-proton collisions 
(Yamazaki et al. 2006). CO observations can detect 
all CO cloud components in the direction of G353.6-0.7, 
therefore help to confirm and distinguish different accel- 
eration mechanisms for emission of TeV 7-rays. E.g. if 
there was a CO cloud overlapping the radio, X and 7- 
ray emissions at the same distance as SNR G353.6-0.7, 
and if a broadened CO emission profile was detected, 
indicating a shock interaction, it would be a strong ar- 
gument to support the SNR interacting with the cloud 
to emit TeV 7-ray emission. However, the existing ^^CO 
observations in the direction of G353.6-0.7 have too low 
resolution (1/8 degree) to give any constraint regarding 
SNR-CO cloud interaction. Further high-resolution CO 
and deep X-ray observations can be used to reveal the 
real radiative mechanism for TeV emission. 

In summary, this new discovery of a radio SNR overlap- 
ping the TeV source HESS J1731-347 provides the best 
case that an old SNR emits TeV 7-rays. 
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